Solitary-wave solutions in binary mixtures of Bose-Einstein condensates under 
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We derive solitary-wave solutions within the mean-field approximation in quasi-one-dimensional 
binary mixtures of Bose-Einstein condensates under periodic boundary conditions, for the case of 
an effective repulsive interatomic interaction. The particular gray-bright solutions that give the 
global energy minima are determined. Their characteristics and the associated dispersion relation 
are derived. In the case of weak coupling, we diagonalize the Hamiltonian analytically to obtain the 
full excitation spectrum of "quantum" solitary-wave solutions. 
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I. INTRODUCTION 

The question of solitary-wave solutions in trapped 
atomic Bose-Einstein condensed gases has received con- 
siderable attention in recent years. Remarkably, solitary 
waves have been created and observed experimentally in 
both one-component [J-Q and two-component systems 
The easiest problem that one can consider is that of 
a quasi-one-dimensional Bose-Einstein condensate which 
extends to infinity with an effective repulsive interaction. 
In this case and within the mean-field approximation, 
the system is described by the nonlinear Schrodinger 
equation, which supports solitary-wave solutions, as was 
shown initially by Tsuzuki Q , and by Zakharov and Sha- 
bat 0. 

If one imposes periodic boundary conditions in a finite 
period of some fixed length, there are significant differ- 
ences. It is no longer possible to work in the thermody- 
namic limit, and it is necessary to impose the particle 
number normalization condition as well as the constraint 
of periodicity. This problem has been addressed within 
the mean-field approximation by Carr et al. Q , for static 
solutions and by the authors of this study for the more 
general case of moving solutions Q . These were shown to 
be Jacobi elliptic functions. The excitation spectrum of 
the many-body problem of a single-component Bose gas 
interacting via a contact potential has also been studied 
by Lieb [10| . This study assumed periodic boundary con- 
ditions and also considered the limit of an infinite period 
Lq and an infinite number of particles TV with N / Lq con- 
stant. The excitation spectrum was found to consist of 
two branches, one of which was later identified as corre- 
sponding to solitary- wave excitation [Tl| . 

The problem of solitary-wave solutions in binary mix- 
tures of Bose-Einstein condensates has also been stud- 
ied theoretically by many people, see, e.g., [l2| - [20j . In 
the case of strictly one-dimensional motion with the two 
condensates extending over an infinite line, the coupled 
system of two non-linear Schrodinger equations that de- 
scribe the two order parameters (within the mean-field 
approximation) has been studied in the context of inte- 



grable systems by Manakov [2l| . 

In the present study we present mean-field solitary- 
wave solutions for a two-component Bose-Einstein con- 
densate with repulsive interatomic interactions for a pe- 
riodic system of finite length Lq. To determine these 
solutions we impose constraints on the particle number 
for each component as well as the constraint of period- 
icity for each species. This formidable set of constraints 
is necessary to determine the solution and the dispersion 
relation, i.e., the energy versus the angular momentum, 
of this system. 

As in the case of a single component, the solitary- wave 
solutions are again found to be Jacobi elliptic functions 
[22l ]. In this case, however, the presence of the second 
component permits several qualitatively different solu- 
tions, i.e., "gray-gray" and "gray-bright" solitons. We 
will argue that the gray-bright solution has the lowest 
energy for particular choices of the phase winding num- 
bers and for interactions which are weak in comparison 
with the kinetic energy of the atoms. The numerical solu- 
tion of the equations which result from these constraints 
also suggests that the choice of the winding numbers does 
not change for moderate couplings. 

This problem is also intimately connected with the 
"yrast" problem, i.e., with the evaluation of the state 
of lowest energy with some fixed angular momentum. 
The yrast state was determined in Refs. }22l.l23| and coin- 
cides [U, [25| with the gray-bright solution found below. 
The connection with the results of Refs. [HHHf, obtained 
directly by minimization of the energy at fixed angular 
momentum, is clarified and gives further support to the 
belief that the yrast states are indeed the gray-bright so- 
lutions. Although gray-gray solutions are also possible, 
the gray-bright solution has a lower energy. The quali- 
tative explanation of this result is simple. The density 
depression found in the gray component is (at least par- 
tially) filled by the bright component. This leads to a 
more uniform total density and thus to a lower energy 
provided that the interatomic interaction is repulsive. 

Finally, in the particular case of weak interatomic in- 
teractions, the complete energy spectrum of the system 
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is determined by diagonalization the many-body Hamil- 
tonian. This is done by identifying certain bilinears of 
the annihilation and creation operators with angular mo- 
mentum operators. Under the assumption of weak cou- 
pling, the interaction terms in the Hamiltonian can be 
re-expressed in terms of Casimir operators and thus per- 
mit the analytic diagonaliztion of the Hamiltonian. 

In the following we first present our model in Sec. II 
and adopt an ansatz that allows us to solve the two cou- 
pled nonlinear Gross-Pitaevskii equations. In Sec. Ill we 
evaluate the energy and the angular momentum of these 
solutions. In Sec. IV we impose the constraints that are 
set by particle normalization and periodicity. In Sec. V 
we consider the nature of these solutions in the limit 
of weak interactions, and in Sec. VI we present numer- 
ical results for our solutions for stronger coupling. In 
Sec. VII we present the results from diagonalization of the 
many-body Hamiltonian. Conclusions and an overview 
are given in Sec. VIII. 
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Making the ansatz 



us = kua + A, 



(0) 



(7) 



(8) 



where k and A are parameters independent of the space 
and time variables, we can integrate these equations. 
With this ansatz, the left side of Eqs. © and ([7]) be- 
come functions of ua and ns respectively that can be 
integrated to yield 
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(lAA + KlAB)n\ 



II. MODEL AND SOLITARY- WAVE 
SOLUTIONS 

Within the mean-field approximation, the order pa- 
rameters of the two distinguishable species (labelled as 
A and B) satisfy the coupled system of the following in- 
tcgrable non-linear equations (Manakov system), 

= -^A + ilAA^A? +1AB\M 2 )1>A (1) 
i^f = -\i'' B + (lBAm 2 +lBB\M 2 )^B, (2) 

where h = 1 and the masses of the two species are as- 
sumed to be equal (and are also set to unity). In addition 
7ij is the matrix element for collisions between species i 
and j. 

The solitary-wave solutions have the form of traveling 
waves with the particle density of each species moving 
with a constant velocity, u, 



vpA = y^I)e l * A(z) e" VAt 

^B = V^e'*^^- 1 " 8 ', 



(3) 
(4) 



where z = x ~ ut, ua and ub are the particle densities, 
and ha and \xb are the chemical potentials of the two 
species. Here x is the spatial variable which is assumed 
to be periodic on the interval < x < Lq. Following 
standard procedures, we separate the real and imaginary 
parts of these equations to find that 



Ca,_ 



riA,B 



(5) 



where Ca and Cb are constants of integration, and also 



1 2 ^_ \C\ 



-2(-u 2 + [i a ^lAB)n\ + E A n A - C\ 



IB A 



1bb)ti 



-2{\u 2 + [i B + -IbaWb + E B n B - C%, 



(9) 



(10) 



where Ea and Eb are integration constants. Consistency 
of the ansatz translates into three equations that relate 
the integration constants and the chemical potentials of 
the two species together with the condition arising from 
the identification of the coefficients of n 3 A 



JAA + KJAB = JAB + K-lBB = 7- 



(11) 



Our ansatz constrains the integration constants, thus re- 
stricting the full solution space of our system of equa- 
tions. In the generic case there are six integration con- 
stants in Eqs. dU and (fTUj) . namely Ea, Eb , ^a, Hb,Ca 
and Cs as well as the two constants n and A arising from 
the ansatz. There are also four consistency conditions 
which reduces the number of free constants to four for any 
given propagation velocity, u. The integration constant 
arising from the integration of Eqs. © and (JTUJ) is not 
included in this counting since it merely corresponds to 
a translation of the solution. However, the solution must 
also satisfy five constraints, namely two constraints of 
particle-number normalization, two phase-matching con- 
straints, and one constraint which sets the period of the 
solution to Lq. In short, there are too many constraints. 

One possible way out of this dilemma would be to 
view the velocity of propagation, u, as a parameter to 
be set by the constraints. This, however, would lead to 
the unphysical result that the velocity of the waves can- 
not be changed without altering the properties of the 
atoms involved. If we ignore the ansatz, u is expected 
to be a free parameter, since we have a total of six in- 
tegration constants (including the chemical potentials) 
and six constraints (i.e., two particle number normaliza- 
tions, two phase matchings and two density matchings.) 
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The restriction on u must be viewed as an artifact of the 
ansatz. 

A more satisfactory way to deal with this problem is to 
fine tune the coupling constants. If the masses of the two 
components are equal, it follows that jab and 7ba are 
trivially equal. If, however, ^ AA = ^ab = Ibb = 7o (i-e., 
if the scattering lengths between the same and the differ- 
ent species are all equal), then the condition of Eq. (fTTj) 
becomes trivial, and we have five free constants and five 
constraints. This allows the velocity of propagation to be 
a free parameter. We thus proceed under the assumption 
that all the coupling constants are equal. Note that in 
this case Eq. (TiTj) implies that 7 = 7 (1 + k). 

It is convenient to factorize the right sides of Eqs. ([9]) 
and (|10|) to obtain the equations 



-72/4 = 27(71.4 - PAi){n A - PA2)(n A - Paz) (12) 
1 27 

-zri!g = — (n B - PBi){n B - p B 2){n B - Pbz), (13) 

A K 

where the roots in these equations are written in ascend- 
ing order. Compatibility between Eqs. ([9]) and (fTO)) and 
Eqs. ([T2[). and (JT3J requires that 



Ca = 1PA1PA2PA3 

n 2 _ 7 

L-B — —PB1PB2PB3- 
K 



(14) 
(15) 



Since both densities have maximum and minimum values 
where 1 their derivatives must vanish, all roots must be 
real. Because of the positivity of n'J, the solution of 
Eq. (|12p is trapped between pai = riA.mhi and pA2 = 
^A,max, i-e., between the minimum and the maximum 
densities of species A. There are two possibilities for 
Eq. (TI"3"T) . For one of these, k > 0, p B \ = "s,min, and 
PB2 = "B.max- For the other, k < 0, p B 2 = n B , m in, 
and PB3 = ^s.max- The first case corresponds to a gray- 
gray solution; the second corresponds to a gray-bright 
solution. 

For the gray-gray solution, we note that Eq. (fT5|) re- 
duces to Eq. $F2§ if 



+ A 
+ A 



(16) 
(17) 
(18) 



The solution of Eqs. (|12[) . and ([13]) can then be expressed 
in terms of Jacobi elliptic functions as 



PB3 = KPA3 + A. 



2K(m)z 

n A = n A ,min + ("A.max ~ n-A,min)sn 2 ( 

2 2K(m)z 

n B = n B . nlin + (n Bimax - n B;min )sn ( 



where 



m 



T^A.max T^A.min f^B .max 7l£i ,: 



PA3 — ™A,min PB3 ~ n B,mh 



m) (19) 
m),(20) 



(21) 



and where the first elliptic integral K{m) satisfies the 
periodicity constraint 



K(m) 



riA- 



0- 



(22) 



Equation (|18|) implies that m is the same in Eqs. (|T9|) and 
(p?D"[) . In this general form, the five independent constants 
in the solution Sp£LC6 clXC 7^ .A, max; HA,mir\.i k, A, and m. 

In the gray-bright case (k < 0), the reduction of 
Eq. ([131) to E q- C2) requires that 



(23) 
(24) 
(25) 



riB,min = ^T^max + A 
KTIA min + A 
PB1 = np A 3 + A. 



The solution of Eqs. (TT2]) and (TJ3]) is now 

tia = riA,min + (nA,m 



n-A, 



, ,2K(m)z , . . 
n )sn 2 ( LJ-| m )(26) 



U 



nB = n B , 



( n B, 



n B ,min)cn ( 



2 2K(m)z 



where 



UAa 



tia, 



n B - 



PA3 - «A,min 



n B, 



PB1 



m),(27) 



(28) 



and K{m) satisfies the periodicity constraint given above. 
Again, the five independent constants in the solution 
space are riA,m& x , n^min, A, and m. 

For both the gray-gray and the gray-bright cases, we 
note that Eqs. (Tl2[) and ([i"3[) have an interesting limiting 
form when k — > — 1 and pA3 — > +00 in such a way that 
(1 + k)pa3 is finite. The fact that pA3 ~ > +00 tells us that 
m — > 0. Hence, the elliptic functions sn(2K(m)z/ Lo\m), 
and cn(2K (m) z / Lo\m) become the regular trigonometric 
functions sin(7rz/L ) and cos(7rz/L )j respectively. The 
periodicity condition of Eq. ([22p assumes the form 



lim 



-1 1 



T 2 

, 
—oloijlA,: 



tia- 



The gray-gray solution then simplifies to 

n A = n A,min + (^A.max — n A,mi 
fl B = n B . m [ n + (n B .nmx — "B.min 

and the gray-bright solution becomes 

1l A = TlA,mm + ("A.max ~ «A,min) sin 2 (-^ 

^0 

n B = n B ,min + (^B.max - ?^B,min) COS 2 ( — ) 

J-'Q 



1 )sm 2 (— ) 

)sm 2 ( — ), 
^0 



(29) 

(30) 
(31) 

(32) 
(33) 



Since the total density is more uniform in gray-bright 
case and since the interaction is repulsive, this solution 
will have a lower energy. The remainder of this paper 
will focus on the analysis of this solution only. 
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III. DISPERSION RELATION: ENERGY 
VERSUS ANGULAR MOMENTUM 

In this section we evaluate the energy and the angular 
momentum of the system in the gray-bright case. We 
begin with the angular momentum, which has the form 

L = ^ / tpA(-i-^)>ipAdx + ^ j i/; B (-i-^)tp B dx 



= 2^/ ( n ^'A + n B^' B )dx. 



(34) 



We can eliminate $^ and $' B using the equations for the 
conservation of particles, Eqs. (JSJ) , to obtain 

^L = u(N A + N B ) + (C A + C B )L , (35) 

where Na and N B are the particle numbers for the two 
species. 

The energy is given by 

/{^lVA| 2 + i|^| 2 + ^7o(|^| 2 + |^| 2 ) 2 }^, 06) 



In the limiting case k — >• — 1 with (1 + k) = 0(m), 
Eq — > 7o(Na + N B ) 2 /(2L ), which is indeed the interac- 
tion energy of a gas of constant density with N A + N B 
particles. In this limit, the remaining energy is 



1 



1 tH 



E-Eq^ -u z N + u(C A + C B )Lo + -j^N. (41) 

Here, we have used the limiting form of the periodicity 
condition, Eq. (f29| , to eliminate quotients of the form 
(1 + n)/m. If the winding number qa = or q B = 0, then 
m = and the velocity becomes u = it/Lq, as it should. 
Since we also have Eq. (p^B]) , it is possible to write 



E-E Q 



2tt2 
T 2 



L 



(Lo/2tt) 



-L. 



(42) 



This is consistent with the equation dE/dL = ft, where 
ft is the angular velocity of the condensate. 



IV. CONSTRAINTS 



which can also be written as 



+ n A (u-i ) 

2 4riA n>A 



(n' B ) 2 



+n B (u-\ ) 2 +"f {n A + n B ) 2 }dx. (37) 

n B 

Using Eqs. (fT2|) -(fT5 |) together with the particle number 
normalization, E can be expressed in the form 

E-E = \u 2 {N A + N B ) + u(C A + C B )L 
7 1 1 

~ o (( 2 W.min + (H W.max)-^ 

2 m m 

7 1 1 

-^-((2 )n B max + (H )n B m[n )N B 

2k m m 



" q (2^ J 4,max^A,mii 



1 



(1 - —)n 2 A in )L 
m 



^(2 7l B,max^S,r 



(l--)n 
m 



2 

B.max 



)L , 



(38) 



where 



Eq 



L /2 



-L /2 



7 2 , 7 2 i 70 / n2i /nn\ 

~ n A + 7T n B + -^\ n A + n B ) }■ (39) 

Z ZK Z 



It is possible to derive a simple form for Eq using the 
normalization constraints (see Appendix 2). In this way 
we get 

Eq = 7o(l + «) 2 n A , min (2N A - n A . mhl L 



2(m + l) , 2(m+l) 



3m 



+7o 



NaNe 
Lq 



3m 

2" Lq 



1 

3m 



3 ^1 2^1 



N 2 



1 iV| 



Lq Lq 2k Lq 



(40) 



The particle number normalization constraints tell us 
that 



L /2 



L /2 



n A ,sdz = Na,i 



(43) 



Using the integrals given in Appendix 2, these equations 
become 

i \^0/i E(m) 

nA,minLo + (llA.max ~ "vi.min) (1 777 s) 

m K (m) 

= N A (44) 

. Lq . E(m) 

ns.min^o + (ns, max - n B min ) — (m - 1 + — -) 

m K (m) 

= iV B ,(45) 

where K(m) and E(m) are the usual elliptic integrals. 
The phase matching constraints, 



r-L /2 

j <&' A B dz — 2irqA.B, 

J-L /2 



(46) 



imply that 



/-La/2 ^ 
dz = 2irq A ,B- (47) 
-L /2 n A,B 

Here q A and (?b are the winding numbers of the two 
species. Equations (|T4")) and (fT5"|) allow us to determine 
CU and Cb as 

C A = 771A, min^A, max "Amin H ( 4 °) 

m 

C B = -nB.minUB,maxPB,max ■ (49) 

K m 
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Carrying out the integrations in Eq. (|47l) and solving for 
the velocity u, we find that 



2wq A 



± 



nA- 



1 



nA.min K(m) 



n(i- 



riA, 



\m) 



.max ^A.min )/' 



(50) 



2^ _ 

^0 V n B,max A'(m) 

X W { r r/K)[n B .max v^B 



n(i 



n B , 



\m) 



n B - 



0/m], (51) 



where n(a|m) is the third elliptic integral. Note that 
the sign ambiguity that appears in Eqs. (|50]) and (fSTj) 
comes from the ambiguity in the sign of the constants 
Ca, Cb- Equations (JSJ) and (fTU)) for the densities ua and 
n b involve only C\ B , hence there is an ambiguity in 
the signs of Ca,b in the solutions. The only place where 
these signs are important are in the phase matching con- 
straints, where the winding numbers also appear. There- 
fore, a solution of the phase matching constraints involves 
not only a choice of the winding numbers qA and qs , but 
also a choice for the signs of Ca,b- The final constraint 
is the periodicity constraint, Eq. (|22[) . encountered ear- 
lier. These five constraints are sufficient to determine 
the solution. 



V. WEAK-COUPLING LIMIT 

In the particular case m — > the periodicity constraint 
of Eq. ([22]) has, to lowest order, a particularly simple form 
(see Appendix 1): 



7T 2 (1 +TO/2) 



riA- 



+ 0(m 2 ). (52) 



Since in general it is not necessarily true that Aha = 
^A,max — fi>A,mm ^ V-^Oj from the above equation follows 



that we should also take the limit n 



T, so that the 



value of lim^^j^L^Q) (1+ k)/to is finite and determines 
An A - 

In this limit, the normalization constraints of Eqs. (|44[) 
and (1451) can be written as 



TlA,: 



n B ,i 



+ «B,min 



An A 

16 
An B 



16 



Nb 
Lq ' 



m + 0{m 2 ) 



(53) 
(54) 



Setting n A ,B = («A,s,max + «A,s,min)/2, we see that 
it-a.b = Na.b/Lq to lowest order in m. 

The phase constants Ca and Cb are given as 



Ca = ±-^y/TlA. T . 



kTIA,] 



1 + m 



n A 

An A 



± — ^A, 

L>Q 



l(1+ 2 m A^ } 



C B = 



IT 

^0 



K n B, 



1 



[1 — — m 



n B 
An B 

n B . 



(55) 



(56) 



Using the asymptotic expansion of n(a|m) and K{m) 
(see Appendix 1), we see that the phase constraints of 
Eqs. (f5T)]) and (f5"Tj) can be written as 



M = 9AT T 2 ) + G(m 2 )(57) 



K ,„ , m \A l B,max'lB,min, o w _„ s 

« = _(2 toT l± T V-^_ ) + 0(m *)( 5 8) 

Here, the signs that appear in the above formulae are to 
be determined by the signs of the phase constants Ca and 
C B (positive, upper sign; negative lower sign), which still 
need to be determined. We also note that these expansion 
formulae are valid only if Aua,b I^a.b ^ Tn since the 
expansions of Eqs. (|55|) and ([56]) are not otherwise valid. 

It is also possible to expand the particle densities near 
m = by making use of the Lambert series of the Jacobi 
elliptic functions (see Appendix 1). This gives us 



n A (z) 



n B (z) 



+ An A [( 

m 



1 

2 ~ 

47TZ 



1 ,2770 

-cos( — )) 

I Lq 



^B,min 

m 



(1 - cos(^))] + 0(m 2 ) 

16 Lq 

a r/1 1 /27TZ.. 

An s [(- + -cos( — )) 



A-7T 7 

_(l-cos(— ))}+0(m 2 ) 

16 Lq 



(59) 



(60) 



The angular momentum given by Eq. ()35l) has the ex- 
panded form 

^-L = u(N A + N B ) ± TTjn A , max nA,min(i + 

±7r^/n B;max n Bimin (l - -m-^-|j51) 
2 Ana 

where it is again necessary to assume that Aua,b 3> m - 
Let us consider now the specific branch (qA,q B ) = 
(0,0). If the velocity u > 0, then Eq. ([4"7|) demands 
that Ca < and C B < 0; hence the lower signs in 
Eqs. ([57)) and ([55]) apply. The only way to realize this 
while still having a common velocity for the two species 
is to have m = exactly. This gives a propagation veloc- 
ity of u = ir(2qA + l)/£o = k/Lq. This means that the 
densities have the form 

. ,1 1 .2to.. . , 
n A {z) = riA.min + An A {- - - cos(—— )) (62) 

Z Z IjQ 

1 1 27TZ 

n B (z) = n B .min + An B (- + -cos(— - )). (63) 

Z Z Lq 
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In this case the angular momentum is 



—L = —(Na + N b )-tt^/1TX; 

J->0 J-'O 



max'*'S,mm* 



(64) 



Recalling that n^max = — ^s.min + N/Lo and that 
«A,min = — TiB,max + N/Lq, it is easy to show that the 
maximum value of the angular momentum in this branch 
is 



1 1 



N 



= -- -VXA-XB- 



(65) 



Making use of the fact that (riB. min + ns !max )/2 = 
Nb/Lq + O(mAns) we note that this maximum is at- 
tained when 7iB,min7iB,max vanishes and this is only pos- 
sible when riB,min = 0. Here, we have written the total 
particle number as N = N a + Nb , the angular momen- 
tum per particle as £ = L/N, and the particle fractions as 
xa,b = Na,b/N. As mentioned earlier, the yrast state 
of a two-component Bosc-Einstein condensate confined 
in a ring trap was evaluated in Ref. [23[ . Based on rather 
general arguments [24|, HH , the present calculation is ex- 
pected to be equivalent to that of the yrast state. Indeed, 
the branch found above corresponds to a portion of the 
first linear branch in the dispersion relation determined 
in Ref. [23| (i.e., for I < xb)- The lowest-energy state 
obtained in Ref. 1231 was found to be 




(co + cie 2 ™ 2 / Lo ) 
(do + die 2 ™^ ), 



(66) 



where 



[x A -t)[l-t) 



£(x B - £) 

x A (l-2£) 

£(x A ~ £) 

x B (l-2£) 



(67) 



for this branch. This gives rise to the following densities 
for the two species 



n B 



These become identical to the densities given by Eqs. (|6"2")l 
and (|6"3"j) if we set 



N A 


N A 


Lo ' 


Lo 


N B 


N B 


L ' 


Lo 



-,d a!i = 



ANa/Lo 



ANb/Lo 



(69) 



This is compatible with the amplitudes of Eqs. (|67|) . 

We consider now the branch (qA, 9s) = (0, 1). It will 
be assumed that the particle wave functions switch con- 
tinuously to this branch as the angular momentum of 



the system increases. When applied to the phase match- 
ing constraints of Eqs. (|4"7} , this continuity demands that 
Ca < and Cb > 0. To 0(m) the phase matching con- 
dition can be satisfied in two ways. One is by having 
m = 0, in which case the densities are exactly as in the 
{IAiQb) = (0,0) case. However, the constant Cb now 
changes sign and becomes positive. This means that the 
angular momentum is given as 



27T 7T 

-f-L = —(N A + N B ) - T-y/nA,m^nA,n 
Lo J^o 



+7r v /n Bimax n s 

.r 



This satisfies the inequality 
1 1 



- - yVxA - XB < £ < X B - 



(70) 



(71) 



This branch corresponds to the remainder of the first 
linear branch in the dispersion relation of Ref. [23| . 
The second way is by having 



y / ^ A ,max^ A ,min y / ^ B ,max^ B ,min 



An A 



An, 



(72) 



to 0(1). This relation together with the ansatz re- 
lations Eqs. (|23[) and (|24[) and the normalization con- 
ditions Eqs. (|53|) and ([54")) to lowest order tell us that 
K = -N B /N A and A = 2N B /L . This is to be under- 
stood as a weak-coupling branch: The periodicity condi- 
tion of Eq. (f52j) demands that An a = 0(m), and this vio- 
lates the condition AnA.B /ua,b ^ fn, unless 70 <C 1/io 
is small. The angular momentum has the form 



Na + Nb^I_V^^a^ 

2 2 Na/Lo tsi \ ) 

Dividing by the total number of particles N, the above 
equation gives 



T l / n A,ma,x n AA 



1 - 2£ 
XA ~ X B 



Na/Lo 

and the propagation velocity becomes 



7T 

To 



2n 



(74) 



(75) 



The minimum and the maximum values of the above an- 
gular momentum are given by the inequality 



x B < I < 



(76) 



This branch can be identified as the first half of the 
curved part of the dispersion relation derived in Ref. [23j j . 
Note that at £ = 1/2 we get u = tt/Lo, which is the ve- 
locity when m is exactly zero. 

The next branch which appears as the angular momen- 
tum increases is given by (qA, Qs) = (1, 0). Here we have 
Ca > and Cb < 0. To order 0(m) the phase matching 
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conditions of Eqs. (|57|) and (|58| may be satisfied in two 
ways. 

One is by again demanding the validity of Eq. (|72j) to 
order 0(1). As before, this leads to k = —Nb/Na and 
to the angular momentum 



L = 



Na + N e 



1 y/TlA 

2 



N A /L 



■(Na-Nb). (77) 



Hence, 



21-1 

- x B ' 



(78) 



and the propagation velocity is given again by Eq. (|75j). 
However, the velocity is now lower than that for m = 0, 
i.e., u = tt/Lq. The minimum and the maximum values 
of the angular momentum are given by the inequality 



1 



<£<x A - 



(79) 



This gives the second half of the curved part of the dis- 
persion relation evaluated in Ref. f23j . 

The other possibility is to set m = 0. In this case the 
angular momentum becomes 



27T 7T 

-T- L = ~r{ N A + N B ) + TTynA.max^A.min 

.min, 

and it satisfies the inequality 



x A <£< 



1 1 



\fx~A 



XB- 



(80) 



(81) 



This reproduces a portion of the second linear branch in 
the dispersion relation evaluated in Ref. [23j . 

Finally the branch (c[At<1b) = (1, 1) appears. Here, we 
necessarily have m = 0, and the angular momentum has 
the form (with Ca > and C B > 0) 



^-L = ^-{N A + N B ) + tt^HZ 



(82) 



In this case the minimum value of the angular momentum 
is 



1 1 



?min = ^ + 2^ XA ~ XB ' 



(83) 



and the maximum value is £ = 1. This describes the 
remainder of the linear part of the dispersion relation 
of Ref. [H, in the interval < I < 1. Beyond t = 1 
the picture repeats itself because of Bloch's theorem (2|| , 
which tells us that an increase of I by an integer can be 
attributed to excitation of the center of mass motion. 



E 0.5 




^ -0.5 



0.5 1 

Ang. mom. per particle / 



0.5 1 

Ang. mom. per particle / 




0.5 
Ang. mom. per particle / 



0.5 1 
Ang. mom. per particle / 




40 



30 



0.5 1 

Ang. mom. per particle / 



20 




0.5 1 

Ang. mom. per particle / 



FIG. 1: The parameters m, n, tla,b, max (solid line), nA,B,mm 
(dashed line), the propagation velocity and the dispersion re- 
lation as functions of the angular momentum per particle I, for 
N A = 160, N B = 40, and j L = 0.2. For 1=1/2 there is a 
node in the density of each species (at different points). Also, 
the velocity of propagation u has a discontinuity at I = xa 
and I = xb- 



VI. NUMERICAL SOLUTION OF THE 
CONSTRAINTS 

It is possible to find numerical solutions to the con- 
straint equations in order to determine k and m. The 
constraints that are crucial in determining these param- 
eters are the phase constraints of Eq. (|T7|) . If we restrict 
ourselves initially to the branch (qa, 9b) = (0, 1), we first 
express u in terms of the angular momentum per parti- 
cle £ and obtain k and m as functions of the angular 
momentum £. This solution is then readily extended to 
the branches {qA,qs) = (0,0), (0, 1) and (1, 1), enabling 
us to plot various observables of the solitary waves as 
functions of £. 

One interesting aspect of the periodic gray-bright solu- 
tion with period equal to Lq, which is expected to be the 
yrast state (i.e., the state of minimum energy for some 
fixed value of the angular momentum), is its size rela- 
tive to Lq. A reasonable measure of this size is the ratio 
of the complex to the real period of the doubly periodic 
Jacobi solution, since the complex period controls the 
exponential decay of the solution. In Fig. 2 we plot the 
ratio of the two periods versus the angular momentum. 
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FIG. 2: The ratio of the imaginary to the real period of the 
gray-bright solution for Na = 160, Nb — 40, and 7o-I/o = 0.2 
as a function of /. The period ratio diverges to infinity as 
£ — > TVs /N — 0.2. When i = 1/2 we have the minimal ratio, 
suggesting that we have the most localized solution. However, 
the ratio is « 0.6 and thus the size of the localized wave is 
still approximately two thirds of the size of the period length 
L . 



This ratio becomes infinite when t — > Nb /N, and it has 
its minimum value when I = 1/2. This suggests that the 
solution is most localized when t = 1/2. However, even 
in this case the period ratio is not close to zero, suggest- 
ing that the yrast state is not very localized but rather 
has a size comparable to Lq even for strong interatomic 
interactions. As seen from Fig. 2 when 70-^0 = 0.2, i.e., 
when the ratio between the interaction energy of the ho- 
mogeneous system and the kinetic energy, NjqLq/(4it 2 ), 
is equal to 1, the minimum value of the ratio of the pe- 
riods is w 0.6. We have also found numerically that the 
corresponding minimal ratio for 70^0 — 2, is w 0.37, 
suggesting that the size of the waves is comparable to Lq 
even for strong coupling. 



VII. DIAGONALIZATION OF THE 
HAMILTONIAN FOR WEAK INTERACTIONS 

For sufficiently weak interactions, 70-^0 << 1> it is 
reasonable to truncate our Hamiltonian to the lowest two 
angular momentum modes only. Doing this, the second- 
quantized Hamiltonian becomes 



, 1 70 . f2 2 

+Aa} ) a\a a 1 -+ 



, t 2 2 iA 2 l2 

+ a[ a 1 + b[ b-y 



2aj6|a 6i + 2a\b\ ) aibQ 
+2a\b\aibi + 2a\bj ) aobi + 2a\ ) b\aibo). 



+2alb f Q a b 



(84) 



We can diagonalize this Hamiltonian by considering the 
algebra of the bilinears of annihilation-creation operators 
appearing in it. We define the operators 

nAo = ajao, "so = b' b 0) "Ai = a\a\, n B i = b\bi, 

J\ = ojai, J B = b\b , J B = blh. (85) 



J A 



a\a . 



The first four operators are the usual number operators. 
The last four can be used to generate two copies of the 
SU(2) algebra. Since [uai — uaq, J a] — 2Ja, it is natural 
to define 

£13 = ^(^1-^0), Jb3 = \(j>>Bi ~ n BQ ). (86) 

These bilinears can be divided into three sets. One set 
consists of the operators HAO + nAi and nso+^si- These 
operators commute with all other bilinears and are hence 
central elements. Their eigenvalues are determined by 
the particle numbers Na and Nb ■ The second set consists 



of the operators Ja3, 
of the operators Jb3, 



J A, and j\. The third set consists 

^SL 

Jb-, and J B . The operators in the 
second set commute with the operators in the third set 
and the operators in each set satisfy the commutation 
relations 



[JA3,B3, Ja.b] = Ja : B, [JA3,B3, J a b\ = — J. 



A,B ' 



and 



[Ja,b, J a b\ = ^<7a. 



3, S3- 



(87) 



(88) 



These are the SU(2) commutation relations, which means 
that the algebra of the bilinears splits into a direct sum 
of two copies of the U(l) algebra and two copies of the 
SU(2) algebra. Note the angular momentum operator 
can be expressed in terms of these operators as 



L = n A \ 



1 



(N a + N b ) + (Ja3 + Jb3)- (89) 



It is now possible to split the Hamiltonian into a central 
part and an SU(2) part, H = He + H , where 



H c = ^{N A + N b ) - ~(N A + N B ) 



2Ln 



Nb) + ^NaNb 

^0 



and 



>a3+Jb3)+-^\JaJ b j <-JaJb—Ja3' 

Lq 



(90) 



~Jb3)- 

(91) 
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The linear part of Hq is can readily be expressed in terms 
of the angular momentum. Since the Hamiltonian is rota- 
tionally symmetric, the quadratic part of Hq must com- 
mute with the angular momentum. This places signifi- 
cant constraints on the form of the quadratic part of the 
Hamiltonian. Indeed, if we set Ja.b = Jai.bi + iJA2,B2 

and J a .b = (Jai,bi,Ja2,B2,Ja3,B3), it is possible to 
rewrite Hq in the form 

2 



2 \L 



(L-N/2)-^-(L-N/2) 2 

~ ~2 x 2 

Jb) 2 ~ J A - J b\- 



70 \ff 



(92) 



N B (N b 



1 ]• 



(97) 

It is now possible to determine the yrast energy. Since 
L = N/2 + ttiab, the value of thab is completely de- 
termined, and the minimum energy is obtained for the 
minimum value of jab- For < L < Nb, —N/2 < 
rriAB < — (Na ~ Nb)/2, and the minimum value of jab 
is |m^s| = N /2 — L. Substituting this value of jab into 
Eq. ([57} yields 



L- 



To 
L 



N(N-l)-L) . (98) 



Since the angular momentum L depends only on Ja3 + 

Jb3, it commutes with J A B and with [J a + Jb) 2 - Also, 

since Ja3,B3 are given in terms of the number operators 
in Eq. (|86|). it is clear that their eigenvalues range from 
—Na,b/2 to Na b/%- This means that we are in the spin 
ja.b = Na,b/% representation of the SU(2) algebra. This 



means that Jab anc i Wa + Jb) 2 are given by 

Nam (Nab 



{Ja 



J A.B 



Jb) 2 



2 V 

Jab{3ab 



I I 



1)1, 



(93) 
(94) 



where / stands for the identity operator and jab ranges 
from \ 3A ~3b\ = (N a - N B )/2 to j A + j B = N/2 as a 
consequence of the usual rules for the addition of angular 
momentum. The eigenvalues of the angular momentum 
operator L are L — N/2 + vhab-, where vhab ranges from 
—jab to jab due to Eq. (|89|) . Therefore, H can be writ- 
ten in the form 



K? 1 



To 



3As{j 



AB 



1) 



Na 
2 



7o_ 
Lo 
Na 
2 



1 



and its eigenvalues are 



E 



7o 



[jAB{jAB + 1) 



2tt 

To 

Na 
2 



Nb 
2 



(L- N/2) - ^-(L- N/2) 
Lo 



Nb 
2 



Na 
2 



Nb 
2 



Nb 
2 



The excited energy levels are given by Eq. (|57| . with 
N/2 - L < jab < N/2. 

For Nb < L < Na the minimum value of jab is by 
(Na — Nb)/2, independent of L. In this case the mini- 
mum energy is 



E 



N B <L<N A 



L 



2n 

To 

~N A N B 



To 
Lo 
N B 



1 



-N(N-l) 



L + NL 



(99) 

are given by 



The energy levels of the excited states 
Eq. (JHTJ with (N A - N B )/2 < j AB < N/2 

For N A < L < N, (N A ~ N B )/2 < m AB < N/2, and 
the minimum value of jAB is tuab = L — N/2. Substitu- 
tion into Eq. (|97|) now gives 



^NaKLkN 



1) - N + L 



(100) 

The excited energy levels are again given by Eq. (|9"T|) with 
L — N/2 < jab < N/2. 



VIII. CONCLUSIONS 



Adding to this the contribution from Hq we get that the 
energy eigenvalues are 



E = -^L 2 
Lo 



^NL- 
Lo 



l^N(N. 



1) 



--^NaN b 
2L 



70 r . ,. 

-r\3AB\jAB 



' 2 V 2 



The issue of finding solitary-wave solutions of the non- 
(?tr)car Schrodinger equation is an old problem with vary- 
ing degrees of difficulty. The most elementary question 
is that of a single component which extends to infinity. 
The case of a single component with periodic boundary 
conditions introduces some interesting complications. In 
the presence of a second component, the similar ques- 
l tions introduce additional complications, as there are 
now two coupled equations. Here we have considered the 
(96) case of solitary-wave solutions in a two-component Bosc- 
Einstein condensed gas, which is confined to a zero width 
(i.e., one dimensional) ring of finite radius, therefore re- 
quiring the imposition of periodic boundary conditions. 

Within the mean-field approximation and with the use 
of a reasonable ansatz for the solution, we have integrated 
the coupled nonlinear equations describing order parame- 
ters to find two analytic solutions which can be expressed 
in terms of Jacobi elliptic functions. 
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This problem is also connected to the determination of 
the yrast state, i.e., the state of lowest energy state solu- 
tion given some fixed value of the expectation value of the 
angular momentum. We have shown explicitly that the 
yrast state for this problem is the gray-bright solution, 
in accordance with general arguments |25j]. The corre- 
sponding phase winding numbers that describe the global 
minima depend on the angular momentum of the sys- 
tem, giving rise in this way to various sectors of the dis- 
persion relation, which nevertheless remains continuous. 
For weak coupling, this is shown analytically, however 
the numerical solution of the constraints suggests that 
the situation does not change qualitatively for stronger 
couplings. 

Going beyond the mean-field approximation, we have 
also diagonalized the many-body Hamiltonian exactly in 
the limit of weak interactions, which allows us to trun- 
cate the Hamiltonian to the two lowest-angular momen- 
tum modes. We have thus managed to derive the entire 
excitation spectrum of this many-body system, which in 
a sense corresponds to a "quantum" solitary-wave solu- 
tion. 

Ideally we would like to find solutions for arbitrary 
masses Ma and Mb, for arbitrary coupling constants 
Iaa, Ibb, and jab, and for arbitrary u. We have found 
analytic solutions by imposing the artificial constraints 
that Ma = Mb and jaa = jbb = Jab- We expect 
that small violations of these constraints would lead to 
new (linear) equations that would be non-singular and 
well-behaved. This suggests that the present constrained 
solutions are broadly representative of all solutions which 
do not violate the constraints violently. 
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where q(m) is the nome function. Its expansion for small 
m is 



. „ m m o, 

aim) = 1 h 0(m A ). 

qy 1 16 32 v ' 



(102) 



Making use of the the expansion of the first elliptic inte- 
gral K(m), 

K(m) = ! (\ + im+^m 2 + 0(m 3 )^ , (103) 
we find that 



,2K(m)z, . . ,7rz. 
sn( ^| m )~sm(— ) + 

m , . ttz ,3nz.. 
16 L L 



(104) 



Similarly, we find 



2K(m)z . 7T z . 

cn( \m) ~ cos( — ) + 

L>0 Lq 

-^(- cos(^) + cosfc) + 0(m 2 ). (105) 

16 Lq Lq 



Expansions similar to Eq. (|103[) also exist for the second 
and the third elliptic integrals, 

E( m ) = I (l - i m - £m 2 + 0(m 3 )\ , (106) 



and 



II(a|m) 



7T 7T771 / 1 



- 1 



32a 



2 VI - a 4a - a 

2 / 6 \ 

+ 6 +C(m 3 ). (107) 
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Appendix 1 



Appendix 2 

We wish to evaluate the integral appearing in Eq. ([39 
In doing this we will use the integrals 



(108) 



K[m) 2 t , v. 2(A'(m) - E(m)) 

srr{u\m)du = — — - — — LL , 

-K{m) 171 

K[m) 4 . , 2[(m + 2)K(m)-2(m + l)E(m)} 

sn (u \m)du = 5 . 

-K(m) 3rn 2 

(109) 



Let us consider the expansion of the gray-bright solu- 
tion when m is close to zero. The Lambert series for the 
Jacobi elliptic functions tell us that, for q(m) close to 1, 



Recalling that the solution ansatz tells us that ub = 
mi a + A, the normalization condition gives A = {Nb — 
kNa)/Lq. Making this substitution in Eq. (f3U)) . we ob- 
tain 



2K{m)z 
sn( \m) 

Lq 



2- 



K(m)\fm 
q{mf/ 2 



; — - sin — ) 



1 



• / 3nz . 
g(ra) J Lq 



,(101) 



Eq = 7o 



(1 + k) 



A* 



L /2 



n A dz 



N A N B 3 N% 



K- 

2 L 



Nr 



1 N% 



Lq Lq 2k Lq 



(110) 
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The particle normalization condition can be used to re- 
duce the remaining integral to 



/Lo/2 i 
sn 2 (: 
-L /2 



" L °/ 2 2l 2 K{m) Z [ L E(jrn) 

— \m)dz = —(1 - —— -) 

Lq m A (m) 



(111) 



3m 2 



(m + 2)-2(m + l)|^ 
A(m) 



3m 



1 + 2(m + 1) 



Na/L 



- "A,: 



(112) 



The normalization condition also enables to eliminate the 
ratio E(m) / K(m) to obtain 



io/2 i 2K(m)z 



sn 4 ( 

-L /2 



m)dz 



Substituting Eqs. ([TTTj) and (|TT2"|) into Eq. (fTTU|) . we ob- 
tain Eq. (|3D]). 
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